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ABSTRACT 


Pressure  and  heat-transfer  data  over  flow-aligned  cylinders  with 
three  nose  shapes  {hemisphere,  flat-face,  and  rounded- shoulder  flat- 
face)  have  been  obtained  at  nominal  Mach  numbers  of  6,  8,  10,  and  19 
and  over  the  Reynolds  number  range  of  0.  009  x  108  to  2. 16  x  1()6, 
based  on  model  diameter.  The  experimental  pressure  and  heat-transfer 
distributions  are  compared  with  theoretical  predictions  and  with 
selected  previously  published  dat,.  from  other  facilities.  The  agreement 
between  the  experimental  results  and  selected  theories  for  the  hemisphere 
cylinder  model  is  good.  However,  the  pressure  and  heat-transfer 
distributions  on  the  other  two  configurations  could  not  be  adequately  pre¬ 
dicted  over  the  entire  model  surface. 
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NOMENCLATURE 


Model  skin  thickness,  ft 

Nose  drag  coefficient  (obtained  by  integrating  the 
experimental  pressure  distributions) 

Pressure  coefficient,  (p  -  p^/q^ 

Maximum  pressure  coefficient,  (ps  - 

Model  skin  specific  heat,  Btu/lbm°R 

Model  diameter,  in. 

Total  enthalpy,  Btu/lbm 

Model  length,  in. 

Mach  number 

Pressure,  psia 

Dynamic  pressure,  psia 

Heat-transfer  rate,  Btu/ft^  sec 

Model  cylinder  radius,  in. 

Free- stream  unit  Reynolds  number,  ft”1  or  in.”1  as  noted 

Reynolds  number  based  or,  cylinder  diameter 

Radius  of  model  cross  section,  in. 

Distance  along  model  surface  measured  from  forward 
stagnation  point,  plus  along  upper  surface,  in* 

Stanton  number,  - - =r-r- 

Pod  U®  (H0  -  Hv/) 
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T  Temperature,  °R 

t  Time,  sec 

u  Velocity,  ft/ sec 

w  Model  skin  specific  density,  j.bm/ft3 

x'  Distance  along  model  axis  measured  from  model  bow  shock,  in. 

x"  Distance  along  model  axis  measured  from  nose- cylinder 

tangency  point,  in. 

A  Bow  shock  stand-off  distance,  in. 

-y  Ratio  of  specific  heats 

6  Angle  between  a  normal  to  the  model  surface  and  the 

free- stream  velocity  vector,  radians 

p  Density,  lbm/ft3 

cr  Angle  between  tangent  to  model  surface  and  free-str  earn 

velocity  vector,  deg 

SUBSCRIPTS 

e  Local  conditions  at  the  edge  of  the  boundary  layer 

o  Stilling  chamber  stagnation  conditions 

SH  Model  shoulder  (tangent  point)  conditions 

s  Model  stagnation  conditions 

w  Model  wall  conditions 

a>  Free-stream  conditions 

SUPERSCRIPT 

*  Critical  conditions  (i,  e. ,  conditions  where  the  velocity  is 

equal  to  the  local  speed  of  sound) 
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The  designers  of  hypersonic  configurations  must  be  able  to  nnake 
a  reasonable  prediction  of  the  flow  field  around  the  body  in  order  to 
estimate  vehicle  performance  and  structural  loadings.  The  use  of  blunt 
.  nosed  configurations  is  dictated  by  the  heating  problems  encountered  at 

hypersonic  speeds.  Therefore,  the  designer  is  faced  with  the  task  of 
predicting  the  flow  field  over  a  blunt  nosed  body  which  is  characterized 
«  by  subsonic,  transonic,  and  supersonic  flow  regions. 

The  solution  of  the  subsonic  and  transonic  blunt  body  problem  has 
been  the  subject  of  considerable  research;  various  blunt  body  solutions 
are  reviewed  by  Van  Etyke  in  Refs.  1  and  2.  The  supersonic  flow  field 
on  a  blunt  body  can  be  obtained  from  the  commonly  used  method  of 
characteristics.  Of  course,  the  validity  of  theoretical  predictions  de¬ 
pends  upon  their  agreement  with  experimental  data.  However,  there  is 
a  definite  scarcity  of  experimental  data  on  blunt  shapes  at  hypersonic 
Mach  numbers  which  can  be  used  to  evaluate  available  theories.  This 
report  presents  data  on  a  cylinder  with  three  types  of  nose  bluntness 
(hemisphere,  flat-face,  and  rounded- shoulder  flat-face)  at  Mach  num¬ 
bers  6,  8,  10,  and  19. 

The  hemisphere  model  is  the  most  commonly  studied  blunt  body  and 
lias  been  designated  as  an  AGARD  calibration  model  (Ref,  3).  Clark  in 
Ref.  4  has  compiled  hemisphere  data  over  a  Mach  number  range  of  i.  8 
to  21  and  developed  empirical  relations  for  the  hemisphere  stagnation 
point  velocity  gradient,  sonic  point,  and  drag  coefficient.  The  flat-face 
model  represents  the  extreme  in  nose  bluntness  and  has  been  investigated 
analytically  by  Gold  and  Holt  in  Ref.  5.  Kemp,  Rose,  and  Detra  (Ref.  6) 
originated  the  rounded- shoulder  flat-face  configuration  in  order  to  pro¬ 
vide  a  more  stringent  test  of  the  similarity-type  solution  used  in  heat- 
transfer  theories. 

The  purpose  of  this  report  is  to  (1)  evaluate  available  theories  and 
discuss  their  range  of  applicability  and  (2)  make  available  new  data  which 
may  be  correlated  with  other  experimental  data. 

SECTION  II 

4  APPARATUS 

2.1  WIND  TUNNELS 

•  *■ 

The  experimental  data  reported  herein  were  obtained  in  five  of  the 
VKF  hypersonic  wind  tunnels;  Gas  Dynamic  Wind  Tunnels,  Hypersonic 
(B),  (C),  (.E),  (H),  and  (F). 
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2.1.1  TunntU  8,  C,  and  E 

Tunnels  B  and  C  are  continuous,  closed-circuit,  variable  density 
wind  tunnels  with  axisymmetrie  contoured  nozzles  and  50-in.  -diam 
test  sections.  Tunnel  B  operates  at  nominal  Mach  numbers  of  6  ana  8 
at  stagnation  pressures  from  20  to  250  psia  and  from  50  to  900  pSia, 
respectively.  Tunnel  C  operates  at  a  nominal  Mach  number  of  10  at 
stagnation  pressures  of  200  to  2000  psia.  Stagnation  temperatures  up 
to  1350°R  in  Tunnel  B  and  1900°R  in  Tunnel  C  are  utilized  to  prevent 
liquefaction  of  the  air  in  the  test  section.  The  above  operating  condi¬ 
tions  result  in  free-stream  unit  Reynolds  numbers  from  0.  30  x  10®  to 
5.  00  x  10®  per  foot  in  Tunnel  B  and  from  0.  30  x  10®  to  2.  35  x  10®  per 
foot  in  Tunnel  C.  Tunnel  C  and  its  associated  equipment  are  shown  in 
Fig.  la  (top)  {Appendix  I).  Details  of  Tunnel  B  are  similar  to  those  of 
Tunnel  C.  The  test  section  tank  and  safety  doors  allow  the  model  to  be 
injected  into  the  tunnels  for  a  test  run  and  then  retracted  for  model 
cooling  or  model  changes  without  interrupting  the  tunnel  flow. 

Tunnel  E,  Fig.  lb,  is  an  intermittent,  variable  density  wind 
tunnel  with  a  flexible-plate-type  nozzle  and  a  12-  by  12-in.  test 
section.  The  tunnel  operates  at  Mach  numbers  from  5  to  8  at  maximum 
stagnation  pressures  from  400  to  1600  psia,  respectively,  and  stagnation 
temperatures  sufficient  to  prevent  liquefaction.  Minimum  stagnation 
pressures  are  about  one-quarter  of  the  maximum  at  each  Mach  number. 

The  maximum  free-stream  unit  Reynolds  numbers  are  15.  6  x  106  and 
6. 5  x  106  at  Mach  numbers  5  and  8.  respectively. 

2.1.2  Tunnels  H  and  F  (Hotshots) 

Tunnels  H  and  F,  Figs,  lc  and  d,  are  conical  nozzle,  hotshot,  wind 
tunnels  with  50-  and  100-in.  -diam  test  sections,  respectively.  These  wind 
tunnels  use  nitrogen  as  the  test  gas  and  operate  over  a  Mach  number  range 
of  18  to  21  and  a  free-stream  unit  Reynolds  number  range  of  0.03  to  0.80x10® 
per  foot.  A  standard  shot  is  made  by  directly  heating  a  small  volume  of 
nitrogen  gas  with  a  high  current  electric  arc.  When  the  electric  arc 
discharges,  the  initially  confined  working  gas  ruptures  a  diaphragm  located 
near  the  throat  of  a  convergent-divergent  nozzle.  The  pressurized  gas 
then  expands  to  the  test  section  and  dump  tank  through  a  conical  nozzle 
(5-deg  half-angle  for  the  50-in.  Tunnel  H  and  4-deg  half-angle  for  the 
ioo-in.  Tunnel  F)  providing  useful  run  times  of  approximately  iiO  to 
100  msec. 

9 

Recent  refinements  in  Tunnel  F  have  advanced  its  usefulness  as  a  test¬ 
ing  device.  The  recent  advancements  in  Tunnel  F  (Ref.  7)  include  a  high 
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spe^d  data  acquisition  system,  a  flow  visualization  system,  high' 

Reynolds  number  testing  in  the  main  test  section  and  at  Mach  15  in  an 
upstream  test  section,  free-flight  testing  techniques,  and  increased 
run  times. 

A  more  complete  description  of  the  VKF  hypersonic  wind  tunnels 
may  be  found  in  Refs.  7  through  12. 

2.2  MODELS  AND  INSTRUMENTATION 

2.2.1  Pro**ur«  Mcdols  and  instrumentation 

The  model  geometry,  the  basic  model  dimensions,  and  a  list  of  the 
facilities  where  the  models  were  tested  are  shown  in  Fig.  2.  Model 
pressures  were  measured  in  Tunnel  B  with  15-psid  transducers,  and  in 
Tunnel  C  with  1-  and  15-psid  transducers,  switched  in  and  out  of  the 
system  automatically  to  allow  measuring  to  the  best  precision.  Model 
pressures  in  Tunnel  E  were  measured  with  5-psid  transducers.  From 
repeat  calibrations,  the  estimated  Tunnel  B  pressure  measurement  pre¬ 
cision  was  ±0. 003  psia.  The  estimated  Tunnel  C  measurement  precision 
was  ±0.  001  psia  for  pressures  less  than  1  psia  and  ±0.  008  psia  for  pres¬ 
sures  greater  than  1  psia.  The  estimated  Tunnel  E  pressure  measure¬ 
ment  precision  was  ±0.  5  percent.  Additional,  information  concerning  the 
instrumentation  systems  of  the  continuous  tunnels  may  be  found  in  Refs.  8 
and  12. 

Because  of  the  short  run  times  in  the  hotshot  tunnels  (30  to  100  msec), 
close-coupled  variable- reluctance  pressure  transducers  have  been  devel¬ 
oped  by  Smotherman  (Ref.  13)  for  the  measurement  of  pitot  and  model 
pressures.  The  estimated  Tunnel  H  and  F  pressure  measurement  precision 
is  as  follows: 

p/ps  0.007  to  0.02  0.02  to  1.0 

Precision  ±10  percent  ±5  percent ' 

2.2.2  Heat-Transftr  Modal*  and  Instrumentation 

The  heat-fransfer  data  n,  +he  continuous  tunnels  (B  and  C)  were 
obtained  by  using  tnin-skin  models  («0.  040  in, )  and  the  transient  technique 
described  in  Ref.  12.  The  thin-skin  models  had  thermocouples  spot- 
welded  on  the  interior  surface,  and  the  wall  thicknesses  were  measured 
at  each  thermocouple  location.  The  hemisphere  and  flat-face  heat- 
transfer  models  were  constructed  of  347  and  310  stainless  steel,  respec¬ 
tively.  Values  of  specific  heats  were  experimentally  determined  from 
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samples  of  the  model  material  for  use  in  the  data  reduction.  The  varia¬ 
tion  of  specific  heats  with  temperature  was  also  taken  into  account  in 
the  data  reduction. 

The  estimated  precision  of  the  heat-transfer  data  in  the  continuous 
tunnels  is  as  follows: 

St/Sts  0.  05  to  0.  2  0.  2  to  1.  0 

Precision  ±15  percent  ±6  percent 

The  heat-transfer  data  in  the  hotshot  tunnels  were  obtained  with 
calorimeter- type  transducers  which  used  thermocouples  as  temperature 
sensors.  Ledford  (Ref.  14)  described  these  transducers  in  detail,  and 
a  complete  description  of  hotshot  instrumentation  and  recording  equip¬ 
ment  was  given  by  Bynum  in  Ref.  15. 

The  estimated  precision  of  the  heat-ti 
tunnels  is  as  follows:, 

St/Sts  0.  009  to  0.  03 

Precision  ±15  percent 

SECTION  ill 
PROCEDURE 


3.1  TEST  CONDITIONS  AND  PROCEDURES 

A  summary  of  the  test  conditions  for  the  present  data,  as  well  as  . 
for  referenced  data,  is  given  in  Table  I. 

Because  tunnel  free-stream  nonuniformitj.es  might  distort  blunt 
model  pressure  distributions  (Refs.  4,  9,  and  16),  the  models  were 
positioned  in  the  continuous  tunnels  at  axial  locations  which  had  rela¬ 
tively  uniform  free-stream  pitot  pressure  distributions.  Also,  in  most 
cases,  the  symmetry  of  the  pressure  and  heat-transfer  distributions 
was  checked  by  data  on  both  the  upper  and  lower  model  surfaces* 

3.2  DATA  REDUCTION 

All  pressure  and  heat-transfer  distributions  presented  are  nondi- 
mensionalized  by  stagnation  point  values  which  are  listed  in  the  data 
tabulations  in  Appendix  II. 


ransfer  data  in  the  hotshot 


0.  03  to  1.  0 
±10  percent 
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In  the  continuous  tunnels,  transient  heat-transfer  data  were  obtained 
by  injecting  the  models  Into  the  airstream  and  recording  model  wall 
temperatures  on  magnetic  tape  at  a  rate  of  20  times  every  second.  A 
digital  computer  was  used  to  fit  a  parabola  through  21  consecutive  tem¬ 
perature  values  centered  at  0  5  sec  after  the  model  reached  the  tunnel 
centerline.  The  temperature-time  derivative,  dTw/dt,  was  then  obtained 
from  the  parabola  and  used  in  the  heating  rate  equation 

q  =  ytbc  dTw/<2t  (1) 

This  equation  neglects  conduction  and  radiation  losses  which  have  been 
estimated  to  be  less  than  1  percent  of  the  convective  heating  rates.  In 
the  hotshot  tunnels,  heating  rates,  q,  were  measured  directly  by  the  use 
of  calorimeter-type  transducers. 

Local  Stanton  numbers  were  obtained  from: 

St  =,  _ 3 -  (2) 

Poo**  <H0-HW) 

In  the  continuous  tunnels,  the  total  enthalpies,  H0  and  Hw,  were  calcu¬ 
lated  using  measured  values  of  T0  and  Tw  and  the  relationship  H  =  cp  T 
which  assumes  cp  =  constant.  The  free-stream  conditions  in  Tunnels 
B  (M,*,  =  6  and  8)  and  E  (Mffl  =  6)  were  calculated  by  assuming  an  isett- 
tropic  expansion  from  the  stilling  chamber  and  the  ideal  gas  relationships 
of  Ref.  17.  The  Tunnel  C  (Mw  =  10)  flow  properties  were  corrected  for 
real  gas  effects  by  using  the  Beattie-Bridgeman  equation  of  state  and  the 
procedures  of  Ref.  18.  In  the  continuous  tunnels,  the  ratio  of  model  wall 
lo  tunnel  stagnation  temperature  (Tw/T0)  ranged  from  approximately  0.  25 
to  0. 50  for  heat-transfer  models  and  from  0.  70  to  1.  0  for  pressure  models. 
The  method  of  determining  flow  conditions  in  the  hotshot  tunnels  is  briefly 
summarized  as  follows:  instantaneous  values  of  pD  and  pg  are  measured 
and  an  instantaneous  value  of  qs  is  inferred  from  a  measurement  of  qw 
using  Lees'  distribution  (Ref.  19)  or  measured  directly  on  a  2-in.  -diam 
hemisphere  cylinder  heal  probe  (see  Refs,  7  and  10).  Velocity,  hence 
enthalpy  (H0),  is  calculated  from  measured  values  of  ps,  qs,  and  the  heat 
probe  radius,  using  Fay-Riddell  theory.  Ref.  20;  Newtonian  pressure 
distribution  near  the  stagnation  point  and  zero  dissociation  are  assumed. 
With  values  of  p0,  p  ,  and  H0  known,  the  remaining  flow  conditions  (M0, 
Re„,  etc. )  are  calculated  as  described  in  Ref.  21.  For  the  short  run 
times  experienced  in  the  hotshot  tunnels,  the  modal  wall  temperature  was 
essentially  constant  at  540°R  (Hw  =  133.  7  Bf.u/lb™);  thus,  the  ratio 
(Tw/T0)  varied  between  0,  075  and  0. 13  yielding  ''cold  wall"  conditions. 
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SECTION  IV 
THEORIES 

This  section  presents  a  very  brief  discussion  of  the  theories  that 
are  compared  with  the  data  in  this  report.  All  theories  used  are  based 
on  nominal  free-stream  Mach  numbers  andy  »  1.4,  unless  otherwise 
noted. 

4.1  MODIFIED  NEWTONIAN  AND  MODIFIED  NEWTONIAN* 

PRANDTL-MEYER 

tn  Ref.  22,  Hayes  and  Probstein  presented  a  thorough  analysis  of 
Newtonian  theoiy  and  developed  the  equation 

Cp  =  2  sin*  tr  (3) 

Lees  (Ref.  23)  suggested  that  for  a  blunt  bod,  with  a  detached  bow  shock 
the  Newtonian  theory  could  be  modeled  to  match  stagnation  point  condi¬ 
tions  by  letting 


Cn  «i  Cp  sin*  a  (4) 

where  t 

Cp  -  SsZSsl 
i'at*  q«. 

Equation  (4)  may  be  written  asi 

=  sin2  er  +  “  cos2  a  (5) 

The  modified  Newtonian  theory  as  used  in  this  report  is  represented  by 
Eq.  (5). 

In  Ref.  24,  Lees  and  Kubota  developed  the  Modified  Newtonian- 
Prandtl-Meyer  theory.  They  showed  that  by  matching  the  Newtonian 
pressure  and  pressure  gradient  with  that  given  by  the  Prandtl-Meyer 
relation,  better  agreement  with  experimental  pressure  data  was  obtained 
in  the  hemisphere  shoulder  region  (1. 0  <  s/R  <  1.  6).  Wagner  presented 
matching  conditions  as  a  function  of  free-stream  Mach  number  in  Ref.  25. 

4.2  NUMERICAL  SOLUTIONS 

Numerical  solutions  using  the  inverse  method  suggested  by  Van 
Dyke  (Ref.  1),  and  matched  with  the  characteristic  solution  at 
Mw  <  1. 05,  have  been  computed  for  the  hemisphere  cylinder  model. 


6 


AEDC«TR*57*T48 


The  inverse  method  starts  with  given  free-sfcream  conditions  and  an 
assumed  shock  shape.  The  flow  equations  from  the  shock  to  the  body 
are  then  numerically  integrated  by  a  marching  technique  until  the 
body  shape  producing  the  assumed  shock  is  determined.  The  shock 
shape  is  iterated  until  the  desired  body  is  produced.  As  pointed  out 
by  Vaglio-Laurin  and  Ferri  (Ref.  26),  the  basic  weakness  of  the  in¬ 
verse  method  seems  to  be  that  slightly  different  shock  shapes,  which 
are  difficult  to  distinguish  from  each  other,  can  lead  to  radically  dif¬ 
ferent  body  shapes.  This  is  one  of  the  reasons  that  numerical  solu¬ 
tions  are  presented  only  for  the  hemisphere  cylinder  model. 

At  high  Mach  numbers  and  sufficiently  low  Reynolds  numbers  (viz. 
hotshot  tunnel  conditions),  boundary-layer  effects  become  very 
important,  as  illustrated  in  Refs.  27,  28,  and  29.  Because  of  the  thick 
boundary  layers  at  these  conditions,  the  "effective”  body  undergoes 
significant  changes  and  can  be  approximated  by  adding  the  boundary- 
layer  displacement  thickness  to  the  original  geometric  body.  In  this 
report,  the  method  given  in  Ref.  27  was  used  to  obtain  the  hemisphere 
cylinder  numerical  solutions  for  the  conditions  in  Tunnels  H  and  F. 

The  inviscid  flow  field  in  an  ideal  (point)  source  flow  was  iterated  with 
a  viscous  boundary-layer  solution  until  there  was  negligible  change  in 
the  "effective"  body  geometry  and  thus  in  the  pressure  distribution  along 
the  model.  When  boundary-layer  effects  were  considered,  the  ideal 
source  flow  field  in  the  inviscid  layer  was  used  for  only  the  character¬ 
istics  solution  and  not  the  Van  Dyke  blunt  body  solution.  When  boundary- 
layer  effects  were  net  considered  (i.  e.  inviscid  predictions),  both  the 
blunt  body  and  characteristic  solutions  were  based  on  parallel  flow. 


4.3  BELOTSERKOVSKII’S  METHOD 

The  Belotserkovskii  method  of  predicting  the  pressure  distribution 
on  a  blunt  body  consists  of  dividing  the  shock  layer,  the  region  between 
the  blunt  nose  and  the  detached  shock  wave,  into  one,  two,  or  more 
strips  depending  on  the  accuracy  desired.  The  equations  of  motion  are 
integrated  from  the  body  to  the  boundary  of  each  strip.  The  boundary 
conditions  are  obtained  from  conditions  on  the  axis  of  symmetry  and  at 
the  sonic  point,  which  is  fixed  at  the  sharp  corner  on  a  flat-face  model. 

The  application  of  Belotserkovskii's  method  to  the  configurations 
investigated  in  this  report  is  limited  to  that  presented  by  Gold  and  Holt 
in  Ref.  5.  They  have  obtained  the  first  approximation  to  the  flat-face 
model  pressure  distribution  at  Mw  -5.8. 
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4.4  VINOKUR'S  SOLUTION 

While  Newtonian  theory  is  remarkably  accurate  for  spherical  noses, 
it  becomes  very  poor  for  blunter  shapes  as  will  be  shown.  Vinokur 
(Ref.  30)  developed  an  approximate  analytic  solution  for  blunt  body  in- 
viscid  hypersonic  flow  based  on  a  constant- density  assumption.  The  use 
of  this  solution  gives  results  which  are  not  expected  to  be  valid  past  the 
sonic  point.  This  method  is  illustrated  in  Ref.  31  where  the  pressure 
distributions  on  a  flat-face  model,  as  well  as  on  other  configurations, 
was  presented  as  a  function  of  the  density  ratio  across  the  detached  shock 
wave.  This  ratio  is  the  only  parameter  which  need  be  specified  when 
using  the  constant- density  solution. 


4.5  BUST  ANALOGY 

.  v 

The  sudden,  concentrated  energy  addition  to  the  flow'  by  a  blunt  nose 
body  at  hypersonic  speeds  may  be  regarded  as  analogous  to  an  explosive 
release  of  energy.  Lukasiewicz  (Ref.  32)  used  this  analogy  and  hyper¬ 
sonic  small  disturbance  theory  to  develop  simple  equations  which  predict 
the  pressure  distribution  and  shock  shape  on  axisymmetric  bodies.  The 
second  approximation  equation  for  the  pressure  distribution  on  axisym¬ 
metric  bodies  which  was  used  in  this  report  is  as  follows: 


0.067 


Moo3\/Cd 

_____ 


+  0.44 


(6) 


The  Cq  values  were  obtained  by  integrating  the  pressure  data  and  the 
x'  values  include  the  shock  stand-off  distance,  A.  Schlieren  or  shadow¬ 
graph  pictures  were  used  to  determine,  the  shock  stand-off  distances, 
and  these  values  are  tabulated  later.  In  most  cases  the  blast  analogy  is 
presented  in  terms  of  the  ratio  p / ps  which  was  obtained  by  multiplying 
Eq.  (6)  by  p^/pg,  values  of  which  are  listed  in  the  figures. 


4.6  LOVE'S  EQUATION 

The  deficiencies  of  blast  wave  theory  in  the  prediction  of  inviscid 
induced  pressures  are  mainly  in  two  regions:  near  the  nose  and  far 
downstream  (viz. ,  p/pa<1.0).  Love  (Ref.  33)  used  the  blast  wave 
pressure  decay  laws  (Ref.  34)  and  attempted  to  correct  these  deficiencies. 
Love's  equation  is: 
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This  expression  assumes  a  known  shoulder  pressure  ratio,  Pgjj/Pg, 
and  decays  from  this  ratio  toward  the  free-stream  static  pressure 
ratio,  pw/ps,  as  the  axial  distance  downetreahi  increases.  Of  course, 
the  question  arises  as  to  what  value  ;of  shoulder  pressure  ratio  should 
be  used.  For  the  hemisphere  model  at  hypersonic  Mach  numbers, 

Love  suggested  a  value  of  Psh/ps  =  9*  ^45  which  ^grees  well  with  most 
of  the  present  data  and  is  the  value  used  in  this  report.  Clark's 
empirical  equation  (Ref.  4)  predicted  shoulder  pressure  ratios  which 
were  4  and  16  percent  higher  than  Love's  value  at  Mach  numbers  19 
and  6,  respectively;  however,  Clark's  fairing  included  viscous  effects 
and  Mach  numbers  considerably  below  6.  '  '  , 

Attempting  to  predict  the  afterbody  pressure  distribution  on  the 
flat-face  cylinder  model  by  Love's  method  points  out  the  main  short¬ 
coming  in  applying  Love's  equation.  That  is,  the  shoulder  pressure 
for  this  model  cannot  be  predicted  and,  in  fact,  the  data  show  that  it 
varies  by  approximately  two  orders  of  magnitude  around  the  sharp 
corner.  For  this  reason.  Love's  equation  is  not  used  to  predict  the 
afterbody  pressures  on  the  flat-face  cylinder  model. 

The  application  of  Love's  expression  to  the  rounded-shoulder  flat- 
face  cylinder  is  also  somewhat  questionable.  However,  since  there  is 
a  well  defined  shoulder  pressure  at  the  tangency  point  of  the  nose  and 
afterbody,  the  data  are  compared  with  Love’s  equation.  From  the 
present  data,  a  shoulder  pressure  ratio  of  Pgjj/ps  =  0.  045  was  used  for 
Mach  numbers  6,  8,  and  10  and  a  value  of  Psh/Ps  =  0-0585  was  used  at 
Mach  number  19  for  the  rounded- shoulder  flat-face  model. 


4.7  LEES' THEORY 


Lees  discussed  laminar  heat-transfer  distributions  over  blunt- 
nosed  bodies  at  hypersonic  speeds  in  Ref.  19  and  developed  the 
equation: 


(t)  fe)  fe) r  _  )  H* 

r  to&vj 4 1  ro-v-sa  v. 

[To  N*'*' \  V  J  ^[\°oo'\uV/,J  J 


(8) 


for  bodies  of  revolution.  Given  a  pressure  distribution,  local  condi¬ 
tions  at  the  edge  of  the  boundary  layer  (denoted  b  *  subscript  e)  can  be 
calculated  by  assuming  an  isentropic  expansion  from  the  stagnation 
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conditions  behind  a  normal  shock  (denoted  by  subscript  s>.  Stagnation 
point  velocity  gradients,  (due/d0)s,  were  experimentally  investigated 
by  Trimmer  in  Ref.  35.  He  correlated  [dMe/d{s/s*)}s  for  blunt 
bodies  at  Mach  numbers  6,  8,  and  10.  In  order  to  use  this  correlation 
parameter,  it  is  necessary  to  modify  the  bracketed  quantity  {}  in  Eq,  (8). 
By  using  perfect  gas  isentropic  flow  relationships  it  can  be  shown  that 


(9) 


For  the  hemisphere  model,  dMe/d  =  dMe/0*d(-^)  where  0*  »  0. 72 
radians.  For  the  flat-face  model,  dMe/d  =  dMe/d  ,  since 
s#  =  R  and  for  the  rounded- shoulder  flat-face  model,  dMe/d^~^  * 
c*(s*)  *  va*ues  dMe/d  used  in  this  report  are  tabulated 


below: 

Model 

dM/d(JL) 

Source 

Hemisphere 

0.97 

Numerical  solution 

Flat-face 

0.29 

Ref.  35;  data 

Rounded- shoulder  flat-face 

0.29 

Ref.  35;  implied  from 
geometry  and  data 
correlation 

SECTION  V 

RESULTS  AND  DISCUSSION 


This  section  compares  the  present  data  with  the  theories  discussed 
in  Section  IV  and  with  previously  published  data  acquired  at  hypersonic 
Mach  numbers*.  The  pressure  data  over  the  models  are  presented  in 
two  groups  of  figures.  The  first  group,  Figs.  3  through  8,  includes  the 


*The  presentation  of  previous  data  is  approximate  since  in  many 
cases  the  data  were  taken  from  figures  without  fine  grids. 
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pressure  distribution  from  the  stagnation  point  to  the  nose-cylinder 
tangency  point.  The  tangency  points  in  terms  of  s/R  for  the  three  models 
are  given  below: 

Model  s/R 


Hemisphere  cylinder  1.57 

Flat-face  cylinder  1.  00 

Rounded- shoulder  flat-face  cylinder  1. 14 

The  second  group.  Figs.  9  through  16,  includes  the  cylindrical  afterbody 
pressure  distributions.  Both  the  nose  and  afterbody  heat-transfer  data 
are  presented  as  a  single  group,  Figs.  17  through  21,  since  there  Is  a 
very  limited  amount  of  afterbody  heat-transfer  data.  The  present  data 

and  the  corresponding  free-stream  conditions  may  be  found  in  Tables  II 
through  VI. 


5.1  MODEL  NOSE  PRESSURE  DISTRIBUTIONS 

Theories  and  data  on  the  hemisphere  model  at  Mach  numbers  6,  8, 

10,  and  19  are  compared  in  Figs.  3a  through  d,  respectively.  Both  the 
numerical  solution  (Van  Dyke  and  characteristics  based  on  inviscid 
parallel  flow)  and  the  modified-Newtonian  Prandtl- Meyer  theories  are  in 
good  agreement  with  the  data.  As  expected  there  was  no  discernible 
Reynolds  number  trend.  In  Fig.  3,  as  well  as  many  of  the  following 
figures,  data  are  presented  from  both  the  upper  and  lower  model  surfaces. 
These  data  are  denoted  by  the  same  symbol  appearing  twice  and  normally 
coincide  as  they  should  for  symmetric  models.  However,  in  some  in¬ 
stances  there  is  a  discernible  difference.  This  difference  can  probably 
be  attributed  to  either  data  precision  limitations  or  free-stream  flow  non¬ 
uniformity  effects.  In  order  to  minimize  these  adverse  effects  and  also 
for  reasons  of  clarity,  data  at  different  Reynolds  numbers  are  averaged 
in  cases  where  there  was  no  significant  Reynolds  number  trend.  An 
average  of  the  present  hemisphere  data  at  each  Mach  number  is  compared 
in  Fig.  4  with  published  data  obtained  after  1957.  A  similar  comparison. 
Fig.  4  in  Ref.  41,  shows  data  obtained  prior  to  1957.  The  numerical 
solution  predicts  a  very  slight  Mach  number  effect  which  is  not  discernible  • 
in  the  scatter  of  the  data. 

Theories  and  data  on  the  flat-face  model  are  compared  in  Fig.  5. 
Modified  Newtonian  theory  does  not  provide  a  reasonable  prediction  of  the 
pressure  distribution  on  this  model  since  it  predicts  a  constant  value  of 
p/pB  over  the  entire  model  face.  In  Ref.  22,  Hayes  and  Probstein  sug¬ 
gested  that  better  agreement  between  the  data  and  Kvwtonian  theory  might 
be  obtained  if  the  shock  angle  were  used  rather  than  body  angle. 
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However,  because  of  the  difficulty  of  predicting  blunt  body  shock  shapes 
and  other  practical  application  problems,  the  body  angle  is  normally 
used  as  was  the  cat*  in  this  report.  Vinokur's  method  yields  the  cor¬ 
rect  trend  for  the  pressure  distribution,  but  shows  an  increasin.  '  mder- 
prediction  in  level  as  the  shoulder  is  approached.  Eelotserkovskii's 
first  approximation  (shown  in  Fig.  5a),  although  somewhat  low  near  the 
shoulder,  gives  the  best  prediction  of  the  pressure  distribution  on  the 
flat-face  model.  Gold  and  Holt  (Ref.  5)  indicated  that  even  better  results 
could  be  attained  with  additional  approximations.  Again,  as  expected,  no 
Reynolds  number  effect  is  indicated. 

i 

The  averaged  present  data  are  compared  in  Fig.  6  with 
Belotserkovskii's  solution  as  well  as  with  previously  published  flat-face 
model  pressure  data.  The  Mffl  <  .10  data  agree  quite  well  with  only  the 
Mach  19  data  indicating  a  few  percent  departure  from  the  data  fairing; 
however,  considering  the  precision  of  the  Mach  19  data  (see  Fig.  5d)  no 
further  conclusions  can  be  drawn. 

A  comparison  of  modified  Newtonian  theory  with  data  obtained  on 
the  rounded- shoulder  flat-face  model  at  Mach  numbers  6,  8,  10,  and  19 
is  shown  in  Figs.  7a  through  d.  Modified  Newtonian  theory  is  the  only 
theory  that  is  readily  applicable  to  this  model;  however,  the  theory  is  as 
much  as  25  percent  high  for  s/R  «  0. 8.  Averaged  data  at  each  Mach  num¬ 
ber  are  compared  in  Fig.  8  for  the  rounded- shoulder  flat-face  model  and, 
as  with  the  other  models,  there  is  no  significant  Mach  number  effect. 

The  lack  of  Mach  number  effects  on  blunt  nose  bodies  has  been  discussed 
by  several  authors  (see  Ref.  44  for  example)  and  is  referred  to  as  the 
independence  principle. 


2.2  MODEL  AFTERBODY  PRESSURE  DISTRIBUTIONS 


Theoretical  and  experimental  pressure  distributions  on  the  hemi¬ 
sphere  cylinder  model  are  compared  in  Fig.  9.  The  method  of  charac¬ 
teristics  shows  the  best  agreement  with  the  data,  although  the  data  are 
generally  a  few  percent  high  which  may  be  caused  by  the  boundary-layer 
growth  along  the  model.  However,  if  boundary-layer  effects  were  the 
cause,  the  pressures  would  decrease  for  a  sufficiently  large  increase  in 
Reynolds  number;  but  this  trend  is  not  detectable  over  the  Reynolds  num¬ 
ber  range  investigated.  Figure  Se  does  shew  that  in  the  shoulder  region 
(s/R  *•  2),  the  viscid  predictions  which  include  source  flow  and  boundary- 
layer  effects  (source  flow  effects  are  negligible  for  the  Tunnel  F  data) 
provide  a  significant  improvement  over  the  inviscid  prediction  shown  in 
Fig.  9d.  In  Ref.  27,  Eaves  and  Lewis  provided  additional  insight  into 
viscous  and  source  flow  effects  on  a  hemisphere  cylinder  at  conditions 
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similar  to  those  of  Figo.  9d  and  e.  The  deficiencies  in  the  blast 
analogy  are  seen  mainly  in  two  regions;  near  the  nose  at  all  Mach  hum-' 
bers  (Figs.  9a  through  d)  and  far  downstream  at  the  lowest  Mach  num¬ 
ber  (Fig.  9a). 

Love’s  equation,  based  on  the  blast  wave  pressure  decay  laws, and 
empirically  matched  at  the  shoulder  and  far  downstream,  is  as  much 
as  25  percent  below  the  data  in  some  regions  (e.  g. ,  Fig,  9c,  s/R  *  8), 

The  sensitivity  of  Love's  equation  to  shoulder  pressure  ratio  (Pgjj/pg) 
is  shown  in  Fig.  9a  where  two  shoulder  pressure  ratios  are  assumed 
(Love's  value  ,>f  0.  045  and  Clark's  value  of  0.  053),  Better  agreement 
is  obtained  with  Clark's  value  ever,  though  it  is  approximately  10  per¬ 
cent  higher  than  the  present  shoulder  data. 

A  comparison  of  an  average  of  the  present  hemisphere  cylinder 
pressure  data  at  each  Mach  number  with  previously  published  data  and 
with  inviscid  characteristics  solutions  is  shown  in  Fig.  10.  The  agree¬ 
ment  anung  the  data  and  with  the  theory  at  similar  Mach  numbers  is 
considered  good.  There  is  an  obvious  Mach  number  effect  on  the  cylin¬ 
drical  afterbody;  however,  the  shoulder  pressures  appear  to  be  inde¬ 
pendent  of  Mach  number  and  do  not  exhibit  the  trend  predicted  by  the 
inviscid  characteristic  solutions.  In  fact,  the  shoulder  pressures  for 
Mm  >  5  are  within  6  percent  of  Love’s  value  of  0.  045.  Also  shown  in 
Fig.  10  are  free- stream  static  pressure?  at  various  Mach  numbers 
calculated  from  an  isentropic  expansion  (7=1. 4).  For  s/R  «  13,  the 
Mach  number  5  and  6  cylinder  pressures  are  slightly  overexpanded,  ‘ 
whereas  the  pressures  at  higher  Mach  numbers  shown  an  increasing 
degree  of  underexpansion. 

*  i 

Comparisons  of  the  blast  analogy  with  flat-face  cylinder  pressure 
data  are  presented  in  Fig.  11.  The  blast  analogy  prediction  is  within 
10  percent  of  the  data  for  s/R  >  5  at  all  Mach  numbers.  For  s/R  <  5, 
the  agreement  is  poor.  Lukasiewicz  (Ref.  32)  points  out  that  the  blast 
analogy  should  not  be  expected  to  apply  in  the  nose  region  since  the 
assumptions  of  hypersonic  small  disturbance  theory  are  violated.  The 
difficulties  of  applying  other  theories  to  this  model  were  discussed4n 
Section  IV.  As  mentioned  previously,  the  shoulder  pressure  (s/R  ^-t)  .. 
varies  by  approximately  two  orders  of  magnitude  because  of  the  extreme 
overexpansion  around  the  sharp  corner,  and  therefore,  the  application 
of  Love's  equation  would  be  somewhat  arbitrary.  There  was  no  de¬ 
tectable  Reynolds  number  effect  on  the  flat-faco  cylinder  for  s/R  >  3. 

The  effect  of  Reynolds  number  on  the  pressure  distribution  for  s/R  <  3 
is  best  depicted  in  Fig.  lie.  Figure  12  shows  a  separation  bubble 
occurring  in  this  region  at  Mw  =  6  and  Re^  =  1.  68  x  10®  in  VKF  Tunnel  E. 
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The  exact  mechanism  of  this  separation  is  not  obvious;  however,  it 
appears  to  be  a  combination  of  leading-edge  separation  (caused  by  the 
sharp  corner)  and  separation  caused  by  the  adverse  pressure  gradient 
provided  by  the  inviscid  flow  field.  These  two  types  of  separation  were 
described  in  Ref.  46. 

Averaged  pressure  data  at  each  Mach  number  are  compared  in 
Fig.  13  for  the  flat-face  cylinder  model.  Data  in  the  region  where  there 
was  a  Reynolds  number  effect  (1.  0  <  s/R  <  2.  3)  were  not  averaged.  In 
general,  the  flat-face  cylinder  afterbody  trends  are  similar  to  those 
exhibited  on  the  hemisphere  cylinder  model  (Fig.  10). 

Comparisons  of  theories  and  rounded- shoulder  flat-face  cylinder 
pressure  data  are  shown  in  Fig.  14.  The  blast  analogy  prediction  is 
within  approximately  10  percent  of  the  data  at  all  Mach  numbers  for 
values  of  s/R  >  5.  Love's  prediction  is  as  much  as  45  percent  below  the 
data  (s/R  «  8),  which  demonstrates  the  magnitude  of  the  erxor  which  can 
be  encountered  by  indiscreet  application  of  Love’s  theory.  Even  though 
there  is  a  well  defined  shoulder  pressure,  it  exists  in  a  region  of  over- 
expanded  pressures  for  6  ^  Mffl  <,  10  and  therefore  the  pressure  increases 
immediately  downstream  while  Love's  equation  predicts  a  pressure  de¬ 
crease.  For  Mffl  =  19  (Fig.  14d),  Love  overpredicts  the  rate  of  pressure 
decrease  immediately  downstream  of  the  tangent  point.  As  with  the  other 
models,  there  4s  no  significant  Reynolds  number  effect  for  s/R  >  3  over 
the  range  investigated.  Averaged  pressure  data  at  each  Mach  number  are 
compared  in  Fig.  15  for  the  rounded-shoulder  flat-face  cylinder  model. 

As  mentioned  previously,  the  cylindrical  afterbody  pressure  distribu¬ 
tions  at  all  Mach  numbers  exhibit  similar  trends.  These  trends  are 
correlated  by  the  blast  analogy  in  Fig.  16.  The  x'/d  coordinate  is  refer¬ 
enced  to  the  bow  shock  location  at  the  stagnation  point  (i.  e. ,  includes  the 
stand-off  distance).  This  choice  of  reference  provided  correlation  nearer 
the  shoulder  than  would  be  possible  if  the  stand-off  distances  were  neg¬ 
lected.  Integrated  pressure  drag  coefficients  and  values  of  shock  stand¬ 
off  distance  ratioed  to  model  diameter  arc  tabulated  in  Fig.  16  for  each 
configuration. 

5.3  HEAT-TRANSFER  DISTRIBUTIONS 

Lees'  distribution  theory  compared  with  hemisphere  cylinder  heat- 
transfer  data  at  Mach  numbers  8,  10,  and  19  is  shown  in  Figs.  17a 
through  c.  For  St/Stg  >  0. 2,  95  percent  of  the  data  were  within  ±6  per¬ 
cent  of  Lees'  theoretical  distribution.  For  St/Sts  <  0. 2,  random  devia¬ 
tions  from  the  theory  of  ±25  percent  were  observed.  Since  the  hemi¬ 
sphere  pressure  data  showed  only  slight  viscous  effects  (primarily  in  the 
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tangency  region),  no  significant  Reynolds  number  effects  would  be 
expected  on  the  heat-transfer  model  and  none  were  noted, 

o- 

-  \  ,A 

Figure  18  compares  an  average  of  the  present  hemisphere  model 
heat-transfer  data  at  each  Mach  number  with  some  previously  published 
data  and  with  Lees'  distribution  at  Mg,  =  8  and  19.  The  agreement 
among  the  data  and  with  the  theory  is  considered  good.  The  scatter  in 
the  data  below  St/Sts  «*  0. 1  is  caused  by  the  difficulty  of  measuring  low 
heating  rates.  Laumarm  discussed  this  problem  in  Ref.  47. 

In  Figs.  19a  and  b,  comparisons  are  presented  between  Lees' 
distribution  theory,  based  on  the  empirical  pressure  data  fairing  of 
Fig.  6  and  flat-face  heat-transfer  data  at  Mach  numbers  10  and  19.  The 
scatter  in  the  data  (s/R  <  1)  for  a  given  Reynolds  number  and  s/R  imply" 
that  the  heat-transfer  distribution  on  this  model  is  more  difficult  to 
measure  than  that  of  the  hemisphere  model  since  this  type  of  scatter  does 
not  appear  in  the  hemisphere  heat-transfer  data.  This  data  scatter  may 
be  caused  by  an  increased  sensitivity  of  the  flat-facq  model  to  free-stream 
flow  nonuniformities  as  compared  to  the  hemisphere  model. 

Forl<s/R<3,  the  pressure  data  on  the  flat-face  cylinder  model 
exhibited  a  Reynolds  number  effect  as  was  shown  in  Fig.  11c.  A  similar 
effect  is  indicated  in  Fig.  19b  which  shows  the  heat-transfer  distribution 
on  the  flat-face  cylinder  model. 

In  Fig.  20  an  average  of  the  present  flat-face  model  heat-transfer 
data  at  each  Mach  number  is  compared  with  the  Leea’  distribution  of 
Fig.  19.  Data  in  the  region  of  the  separation  bubble  were  not  averaged. 
For  s/R  <  1,  93  percent  of  the  data  are  within  +8  percent  and  -1  percent 
of  Lees'  predictions. 

Lees'  distribution  theory  is  compared  with  rounded-shoulder  flat- 
face  cylinder  heat-transfer  data  at  M„  =  19  in  Fig.  21.  Lees'  distribu¬ 
tion  Is  based  on  the  pressure  data  fairing  of  Fig.  !4d  and  is  approxi¬ 
mately  30  percent  above  the  data  for  s/R  >  1.  Kemp,  Rose,  and  Detra 
(Ref.  6)  implied  that  their  theory  provides  a  better  prediction  of  the  heat- 
transfer  distribution  on  this  model.  Unfortunately,,  detailed  calculations 
of  their  theory  for  conditions  corresponding  to  the  present  data  ate  riot 
available. 

SECTION  VI 

CONCLUDING  REMARKS 

Pressure  and  heat-transfer  data  over  flow- aligned  cylinders  with 
three  nose  shapes  (hemispherical,  flat-face,  and  rounded- shoulder 
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flat-face)  have  been  obtained  at  Mach  numbers  of  6,  8,  10,  and  19  and 
over  the  Reynolds  number  range  of  0,  0i)9  x  10®  to  2. 16  x  10®,  based  on 
model  diameter.  The  experimental  pressure  and  heat-transfer  distribu¬ 
tions  were  compared  with  previously  published  data  from  other  facilities 
and  with  applicable  theories.  The  agreement  among  the  data  from  the 
various  facilities  is  considered  very  good.  The  limitations  and  pdnge  of 
applicability  of  the  theories  discussed  showed  that  considerable  care 
must  be  exercised  when  attempting  to  predict  the  pressure  and/or  heat- 
transfer  distributions  on  even  such  basic  configurations  as  those  in¬ 
vestigated  in  this  report.  The  theories  investigated  were: 

1.  Modified  Newtonian 

2.  Modified  Newtonian  Prandtl- Meyer 

3.  Van  Dyke  and  characteristic  solution  (inviScid)  - 

4.  Van  Dyke  and  characteristic  solution  (viscid) 

5.  Belotserkovskii's  Method  ' 

6.  Vinokur's  Solution 

7.  Blast  Analogy 

'  8.  Love's  Equation 

9.  Lees'  theory  based  on  a  theoretical  pressure  distribution  " 

10.  Lees'  theory  based  on  an  empirical  pressure  distribution 

A  summary  of  the  comparisons  between  the  present  data  and  the  above 
theories  is  given  in  the  following  table. 


Model  Regions 

'J 

Configurationsj 

Nose 

0<s/R^0.  9 

Shoulder 

6.9  <s/R  <5 

Afterbody 

5  <  s/R  <5 15 

Pressure 

Heat 

Transfer 

Pressure 

Heat 

Transfer 

Pressure 

Heat 

Transfer 

Hemisphere  Cylinder 

/Sv 

\&  va> 

/a\ 

'VS 

w 

Yo\ 

>as 

TD  ; 

Flat-Face  Cylinder 

® 

® 

NA 

NA 

ID 

Rounded-Shoulder 
Flat-Face  Cylinder 

N 

N 

N 

N 

<D 

ND 

16 


A8tJC*TS*67»U# 


Legend;  <x)  Indicates  that  theoretical  prediction  number  s  (from 
above)  was  within  ±10  percent  of  the  VKF  data  at  all 
Mach  numbers  over  the  entire  s/R  range  listed. 

N  Indicates  that  none  of  the  theories  investigated  were 
within  ±10  percent  of  the  VKF  data.  * 

NA  Indicates  that  none  of  the  theories  investigated  were 
applicable.  " 

ND  Indicates  no  data 
ID  Indicates  insufficient  data 
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